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A generating function for the A/-soliton solutions of the Kadomtsev-Petviashvili II equation 

Sarbarish Chakravarty and Yuji Kodama 

ABSTRACT. This work describes a classification of the yV-soliton solutions of the Kadomtsev-Petviashvili II equation in 
terms of chord diagrams of N chords joining pairs of 2N points. The different classes of A'-solitons are enumerated by the 
distribution of crossings of the chords. The generating function of the chord diagrams is expressed as a continued fraction, 
special cases of which are moment generating functions for certain kinds of ^-orthogonal polynomials. 
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1. Introduction 

The classical theory of enumerative combinatorics has indeed a far-reaching scope, encompassing disparate areas 
in mathematical, physical, biological and social sciences. Combinatorial entities such as permutations, partitions, trees, 
lattice paths, graphs and their various enumerations find applications ranging from econometrics, DNA structures, and 
statistical mechanics to coding theory, knots and enumerative algebraic geometry. The purpose of the present note is 
to elaborate on a somewhat unexpected relationship between a classical combinatorial problem studied by Touchard in 
the 1950s and the classification of a special class of solitary wave solutions (solitons) of an exactly solvable nonlinear 
partial differential equation discovered some 20 years later. This nonlinear wave equation, known after its discoverers 
as the Kadomtsev-Petviashvili (KP) equation, 



(1.1) 



3 / ^du 
dx \ dt 



d 3 u 



du 

DM — 

ox 



d 2 u 



describes the evolution of small-amplitude, weakly two-dimensional solitary waves in a weakly dispersive medium 
|18|. Depending on the sign of a, there are two versions of the KP equation namely, KPI and KPII. Throughout 
this article we consider Eq. dl.U with a = +1, which is the KPII equation. The function u = u(x,y,t) is the rescaled 
amplitude of the wave-form. The KP equation arises in many physical settings including water waves and plasmas 
(see e.g. |[T5l for a review). It is a completely integrable system with remarkably rich mathematical structure which is 
well-documented in several monographs lfll H4ll22l . 

In this article, we consider a family of real, non-singular solitary wave solutions of the KPII equation, known as 
the Af-soliton solutions. At any given time t, these wave-forms are localized along certain lines in the xy-plane, and 
decay exponentially everywhere else. In the generic case, they form a pattern of N intersecting straight lines as \y \ — > °° 
in the xy-plane, whereas in the near-field region the N lines interact to form intermediate lines and web-like structures 
as shown in Fig. II. II The simplest kind of such solution is the 1-soliton, which is a constant amplitude wave localized 
along a line in the xy-plane, and traveling with uniform velocity perpendicular to the line. For N > I, the asymptotic 
form of the iV-soliton solution coincides with N 1 -solitons along different directions, as |y| — > °° and uniformly in t . 
For this reason, these solutions are often referred to as the line-solitons. 

Several researchers lPP2lH3ll23ll24ll3H as well as the authors JSjEiBSlSl 6 1 have studied the soli ton solutions of 
KPII. The general line-soli ton configurations called the (N- ,N+) -solitons consist of N- line solitons as y — > — °o and 
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N+ line solitons as y — > °° [ 3, 6 1 . The Af-solitons correspond to the special case when A^_ = N+ = N, and when the N line 
solitons as y — > °° are pair- wise identical (as wave-forms) to the N line solitons as y — » — °°. An interesting feature of the 
Af-soliton solutions is the fact that these solutions can be essentially (i.e., up to space-time translations) reconstructed 
from the asymptotic data alone, comprising N pairs of amplitudes and directions associated with the N line solitons 
as \y\ — > °° HJ[6j[7). As a direct consequence of this result, it is possible to classify all the A^-soliton solutions into 
(2N — 1)!! distinct equivalence classes, corresponding to the ways of partitioning the integer set {1,2, .. . ,2N} into 
N distinct pairs. The purpose of this paper is to extend our studies to a characterization of the A^-soliton solutions 
according to soliton interaction patterns, and give a classification of such interactions in terms of certain partitions 
(perfect matchings) of the integer set { 1 , 2, . . . , 2N}. We highlight some interesting connections between A^-solitons of 
KPII on one hand, and combinatorics of chord diagrams and g-Hermite polynomials on the other. 




Figure 1.1: N-soliton solutions of the KPII equation illustrating different spatial interaction patterns: (a) 2-soliton solution, (b) 
resonant 3-soliton solution, (c) partially resonant 4-soliton solution. Here and in all following figures, the horizontal axis is x, 
vertical axis is y, and the graphs show contour lines of \nu(x,y,t) for fixed /. 



2. Background 

In this section we give a brief overview of the chord diagrams with 2N points, as well as the line-soliton solutions 
of the KPII equation. The aim is to underscore the connection between these two seemingly disjoint mathematical 
objects. In particular, we illustrate how a KPII line-soliton can be represented by a set of index pairs, which leads 
naturally to the construction of a chord diagram on an integer set. We shall use this construction in Section [3] to 
identify each line soliton of an A^-soliton solution by using a chord joining two specific points among 2N points, and 
study the soliton interaction patterns in terms of such chord diagrams. 

2.1. Chord diagrams. Let us first describe a chord diagram consisting of N chords. Consider a partition of the 
integer set [2N] := {1,2,..., 2N} into N distinct 2-element blocks or pairings 

Vn'=%,in}, l<in<j„<2N, n = l,2,...,N, 

such that [2N] is a union of the blocks pi,P2, • • • , Pw- In combinatorics, such a partition is referred to as a (perfect) 
matching of [2N]. We will denote the set of all matchings of [2N] by (M^. The total number of matchings in H 'n is 
given by 

|<MV| = l-3-5...-(22V-l)=: (22V- 1)!!. 
A standard way to represent a matching X of Mn, is to mark 2N points on a line from left to right labeled by 
1,2,... ,22V, and join the two points of each pairing p by a semicircular arc above the line. The resulting diagram 
(see e.g. Fig. 12. U is called a linear chord diagram, whereas a chord diagram would correspond to labeling the 2N 
points in a clockwise manner on a circle, and joining the two points of each pairing by a chord. 

Without loss of generality, the smallest integer i„ from each pairing of X € !Mn can be arranged in a strictly 
increasing order 1 = i\ < i% < . . . ijy < 2Af — 1. However the j n 's are not ordered in general. 

Definition 2. 1 . Let p r and p v with r < s be distinct pairings (equivalently, a pair of chords). Then, 
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1 2 3 4 5 6 7 8 

Figure 2. 1 : A linear Chord diagram 



(a) p r and p s form an alignment or a O-type configuration if i r < j r < i s < j s . That is, the pairs do not overlap. 

(b) p r and p s form a crossing or a T-type configuration if i r < i s < j r < j s . That is, the pairs partially overlap. 

(c) p r and p s form a nesting or a P-type configuration if i r < i s < j s < j r . That is the pairs completely overlap. 



The meanings of O-, T-, and P-type of configurations will be clear later in Section [3] when we discuss the N- 
solitons of KPII. In Fig. 12. 11 the pairing [4,6] forms an alignment (O-type configuration) with [1,3]; a crossing with 
[5,8]; and a nesting with the pairing [2,7]. Furthermore, the total number of (pairwise) crossings in Fig. 12. II is 3; 
they occur between the pairs [1,3] and [2,7]; [2,7] and [5,8]; and [4,6] and [5,8]. Similarly, there is one nesting, 
{ [2, 7] , [4, 6] }, and two alignments, { [1 , 3] , [4, 6] } and { [1 , 3] , [5 , 8] }. 

It should be clear that the of alignments, crossings and nestings for any X 6 Mn must add up to the total number 
of pairwise chord configurations, i.e., N(N — l)/2. One of the earliest results ifTOl [341 in the enumeration of chord 
diagrams is that the number of matchings in Mn with no crossings is given by the jV th Catalan number Cn — ( 2 ^) , 
which appears in many combinatorial problems (see, e.g., 11321 ). Similarly, the number of diagrams in 94n with no 
nestings is also given by Cn- The problem of counting the elements of 5Ww according to the number of pairwise 
crossings of chords was considered by Touchard [33 1, who gave an implicit formula for the enumerating generating 
function in terms of continued fractions. Subsequently, Riordan [28 1 derived a remarkable explicit formula for the 
generating function based on Touchard's work. If cr(X) denotes the number of crossings of the elements G Mn, then 
the generating function by the number of crossings is defined via the polynomial 

Mi) ■= E i cr(x) > ^ cr ( x ) ^ j n ( n - ! ) > 

in the variable q with positive integer coefficients. The Touchard-Riordan formula for F^{q) is 



1 v 



(2.D F N {q) = ^——^ £(-!)" 



2N \ ( 2N 
N-n) ~ \N-n-l 



The first few polynomials are 

Fi(q) = \, F 2 (q)=q + 2, F 3 (q) = q 3 + 3q 2 + 6q + 5 , 
and it easily follows from Eq. (12. U that the number of non-crossing diagrams is given by 

{2N\ ( 2N \ 1 f2N^ 



Nj \N-lJ N+1\N 

which is the Catalan number Cn mentioned earlier. However, the Touchard-Riordan formula is somewhat mysterious, 
in that it is not obvious from Eq. (12. Il l that FN(q) is in fact a polynomial in q of degree A^(A^ — l)/2 as implied by 
its combinatorial origin, or that Fn(1) = \Mn\ = (2N — 1)!!. These assertions follow only after detailed analysis of 
Eq. ( 12. U [ 28 1 (see also OH ). A purely combinatorial proof of the Touchard-Riordan formula also appeared in Ref. [ 25 1 
(see also [19]), and its relation to g-Hermite polynomials was investigated in ref. |H71 . 



2.2. The T-function and line solitons of KPII. It is well known (see e.g. J30j[l4|) that the solution u(x,y,t) of 
the KPII equation is given in terms of the T-function z(x,y,t ) as 

a 2 

(2.2) u(x,y,t) = 2-^j]nz(x,y,t). 
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We consider the class of solutions whose x-functions are given by the Wronskian determinant form, i.e., 



(2.3) T(x,y,0=Wr(/i,...,/ JV )=det 



( h h ■■■ h 
f'i fi ••• In 



\ AN-l) ,(JV-1) ,(AT-1) 
Vl J 2 Jn 



with fn = d-ifn/dx-i, and where the functions {f n }^ = \ form a set of linearly independent solutions of the linear system 

d ± = d lL d l = d 2l 

dy dx 2 ' dt djc 3 

The soliton solutions of KPII can be constructed from Eq. (12.3b by choosing a finite dimensional solution for each 
function /„ (x,y, t), namely, 

M 

(2.4) fn(x,y,t) = £ a nm e Qm , n = 1,2, . . . ,AT, 

m=l 

where Q m (x,y,t) — k m x + k^y + kf n t + 0^, , m = 1,...,M, are phases with real distinct parameters k\,k 2 . . .kM, and real 
constants B®,... ,Q M . The constant coefficients a nm define the N x M coefficient matrix A :— (a nm ). The simplest 
example is the 1-soliton solution with N = 1 ,M ~ 2, and T = f\ = e 6 ' + e® 2 , for which 

d 2 1 1 

u(x,y,t)=2j 1 lnz = -{h - k 2 ) 2 sech 2 -{Q x - 9 2 ) . 

This 1-soliton solution describes a plane traveling wave-form with constant amplitude (k\ — kj) 2 /2. For fixed t, the 
wave-form is localized in the (x,y) -plane along the line L : Q\ — 02 whose normal has the slope c = k\ +k 2 . The 
solution is characterized by two physical parameters, namely, the soliton amplitude parameter a = \k\ — &2I and the 
soliton direction parameter c = k\ + k 2 - 

In the general case, substitution of Eq. ( 12.4-b into the Wronskian of Eq. ( I2.3K and subsequent development of the 
resulting determinant via Binet-Cauchy formula, yields the following explicit form of the T-function: 

(2.5) x(x,y,t) = A(mi,...,m N ) exp[Q(mi,...,m N )] Yl (km r -k ms ), 

l<mi<".<mjv<M l<s<r<N 

where A (mi, . . . ,111^) is the N x N maximal minor of A obtained from the columns 1 < m\ < ■ ■ ■ < < M, and 
8 (mi, . . . , : = + • • • + Qm N is a phase combination of (out of M) distinct phases. Note that the transformation, 
G : A — » A 1 := GA, G G GL(A r ,K), amounts to an overall rescaling of the minors A (mi,. . . ,m^), and hence, of the 
T-function in Eq. (12.5b : i.e., x — > x' = det(G)x. Since such a rescaling leaves the solution u(x,y,t) in Eq. ( 12.2b invari- 
ant, it is possible to reduce the coefficient matrix A to reduced row-echelon form (RREF) by Gaussian elimination. 
Throughout the rest of this article, the coefficient matrix A will be assumed to be in RREF. 

The solutions u(x,y,t) resulting from the x-function in Eq. (12. 5b are singular for arbitrary choices of the parameters 
{k n W = \ and the matrix A. To avoid such singularities, which correspond to the zero-locus of the x-function, one needs 
to impose certain positivity conditions. 

Condition 2.2 (Positive definiteness of x). 

(a) The phase parameters are distinct, and are ordered as k\ < k 2 < ... < kM- 

(b) The N xM coefficient matrix A satisfies rank(A) = N, and M > N. 

(c) All non-zero N xN minors of A are positive. 



Remark 2.3 The matrices satisfying Condition 12. 21 c) above, are called totally non-negative (TNN) matrices. The 
classification of the (N-,N+) -soli ton solutions is thus given by the classification of the N x M TNN matrices A in 
RREF. From a more geometric perspective, each TNN matrix parametrizes a unique cell in the TNN Grassmannian 
Gr + (N,M) (see e.g. Il26l ). and the classification of the soliton solutions corresponds to a further refinement of the 
Schubert decomposition of Gr(N,M) into TNN Grassmann cells (see [20 1 for the case M — 2N). The refinement 
is given by a classification of the coefficient matrix A whose N x N minors A(m\ ,.. . ,m#) represent the Pliicker 
coordinates of Gr(N,M). It should be noted that each (A r _,A r +)-soliton solution corresponding to a TNN matrix can 



be parametrized by a chord diagram [6, 8|. The geometric structure of this classification will be discussed in a future 
communication Q. 

In Condition 12.21 the distinctness assumption on the set of phase parameters ensures that the set {<? e '"}^ =1 is 
linearly independent, while rank(A) = N implies that the set of functions {/ n }„ =] is linearly independent. Also when 
M = N, the sum in Eq. ( 12.5b reduces to a single exponential term with a phase combination that is linear in x. The 
logarithm of such a T-function is annihilated by the second derivative in Eq. (12.2b , leading to the trivial solution 
u(x,y,t) = 0. Thus the condition N < M guarantees non-trivial solutions of KPII. Condition 12. 2f c) together with the 
ordering k\ < k% < . . . < Icm make the sum in Eq. ( 12.5b totally positive. As a result, x(x,y,t) is a positive function 
on R 3 , and the resulting solution u(x,y,t) of KPII is non-singular, bounded and positive definite. Furthermore, the 
asymptotic analysis of the T-function in Eq. d2.51 > reveals that for any given value of t , there exist a set of lines given 
by {Lij : 0, = 8/, i < j} in the (x,y)-plane, such that 

(2.6) u(x,y,t) ~ \{kj -£,) 2 sech 2 ±(8; - e< + 8y) , 

along each Ly either as y — > °°, or as y — > — oo. Equation ( 12.61 ). which has the same form as the 1-soliton solution 
defines an asymptotic line soliton along Ly associated with the solution u(x,y,t). Each line soliton which is parallel 
to the line L ;/ - has the parameters ay = \ki — kj\ for the amplitude and en = ki + kj for the direction normal to the line 
Lij. Hence, we denote each line soliton by the index pairing p = labeling the line Ly. Note from Eq. (12.41 ) that 
the indices i,j labeling the phases 0, , 7 in Eq. ( 12.61 ) also label a pair of distinct columns of the coefficient matrix A. 
Due to this connection, it turns out that the pairing: p = [i, j], 1 < i < j < M, of the line solitons can be determined 
from the structure of the coefficient matrix A which is in RREF, and satisfies the following irreducibility condition. 

Condition 2.4 (Irreducibility) 

(a) Each column of A contains at least one nonzero element. 

(b) Each row of A contains at least one nonzero element in addition to the pivot. 

Recall that, for anffxM matrix in RREF, the leftmost non-vanishing entry in each nonzero row is called a pivot, 
which is normalized to unity. The index pairs [i, j] of the asymptotic line solitons shown in Eq. ( 12.61 ) are then given by 
the following technical result, which is proved in Ref. J5). 

PROPOSITION 2.5. Let the sub-matrices X[ij] andY[ij] of A be defined in terms of their column indices as 

X[ij] : 1.2 i 1../ • 1 V/ Y i j :=[i+l,. ..j - 1] . 

Then, necessary and sufficient conditions for an index pair [i,j] to specify an asymptotic line soliton, as in Eq. ( 12.61 ), 
are the following rank conditions. 

(i) Each line soliton as y — > °° is labeled by a unique index pair [e„,j„] with e n < j n , where {e n }„ =1 label the 
pivot columns of A. Moreover, ifrank(X[e n j n ]) =: r n , then r„ < N — 1 and 

rank(X[e„;„]|e„) = Taak(X[e n j n ]\j„) = mnk(X[e„j„}\e„J n ) = r„ + l. 

(ii) Each line soliton as y — > — °° is labeled by a unique index pair [i n ,gn] with i n < g n , where {gn}n=i label 
the non-pivot columns of A. Moreover, //rank(F [i n g„] ) —'■ s„, then s„ <N — 1 and 

Twk(Y[i n g n ]\i n ) = rank(Y[i„g„]\g„) = rank(Y[i„g n ]\i„,g n ) = s n + l. 

Here (Z\m, n) denotes the sub-matrix Z of A augmented by the columns m and n of A. 

It should be clear from the above result that the (A^_ , ) -soliton solution of KPII generated from the T-function in 
Eq. (12.5b has exactly = asymptotic line-solitons as y — » °° and = M —N asymptotic line-solitons as y — » — oo. 
From the T-function data consisting of M distinct phase parameters k\ , . . . ,&m and a matrix A satisfying Condition l2.4l 
Proposition 12 . 5 1 provides an explicit way to identify all the asymptotic line solitons of the corresponding solution of 
the KPII equation. We illustrate this method with an example. 

Example 2.6 Consider the solution u(x,y, t) generated by the T-function of Eq. (12.3b in the case N = 2 and M — 4, with 
4 real parameters k\ < ki < k^ < k^, and 

f2=e°i+e°i, A=(l J J 
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The pivot columns of A are labeled by the indices \e\,e-i\ = {1,2}, and the non-pivot columns by the indices 
{8I182} = {3,4}. According to Proposition 12.51 the number of asymptotic line solitons is N+ = N- = 2. They 
are identified by the index pairs [1 , 71] , [2, 72] as y — > 00, for some 71 > 1 and 72 > 2; and by the index pairs [i\ , 3] , [12 , 4] 
as y — > — °°, for some i\ < 3 and i% < 4. We first determine the asymptotic line-solitons as y — > °° using the rank 
conditions prescribed in Proposition 12. 5 ID - For the first pivot column e\ = 1; starting from j = 2 and then repeat- 
edly incrementing the value of 7 by unity, we check the rank of each sub-matrix X[l j]. Proceeding in this way, we 
find that the rank conditions are satisfied only when 7 = 4: X[14] = 0. So, rank (X[ 14]) = < N — I. Moreover, 
rank(X[14]|l) — rank(X[14]|4) = rank(X[14]|l,4) = 1 since columns 1 and 4 are parallel. Thus, the first asymptotic 
line soli ton as y — » °° is identified by the index pair [1,4]. For e^ — 2, proceeding in a similar manner we find that 
7 = 3 does satisfy the rank conditions since X [23] = (q ~} ) is of rank 1 =N— 1, and rank(X[23] |2) = rank(X[23] |3) = 
rank(X[23]|2,3) = 2. Therefore, the asymptotic line solitons as y — > °° are identified with the index pairs [1,4] and 
[2,3]. 

We next consider the asymptotics for y — ► — 00. Starting with the non-pivot column gi = 3, we apply the rank 
conditions in Proposition ^. 5f ii) to the column ; = 2. Then, we have Y[23] = 0, and rank(T [23] |2) = rank(T[23]|3) = 
rank(F [23] |2, 3) = 1 . Hence, the pair [2,3] identifies an asymptotic line-soliton as y — » — 00. For g2 = 4, we consider 
i = 1,2,3 and find that the rank conditions are satisfied only for i = 1. In this case, Y[14] = ( j j), so rank(T[14]) = 
1 = N- 1 andrank(T[14]|l) =rank(7[14]|4) =rank(T[14]|l,4) = 2. Thus, the index pair [1,4] identifies the other 
asymptotic line-soliton as y — » — 00. In summary, both pairs of asymptotic line solitons as y — > ±°° are labeled by the 
index pairs [1,4] and [2,3]. This is an example ofaP-type 2-soliton solution (see Section|3]l, and is shown in Fig. 13. II 



It should be emphasized that in general N- ^ N+, and that even in the case jV_ = N + , the line solitons as y — > °° 
are in general distinct from the line solitons as y — > — °o in both amplitude and direction. In this article, we restrict our 
discussions primarily to the A^-soliton subclass of the line-soliton solutions of KPII. 

Definition 2.7. Let S + := {[e„,j„]}^ =l and S_ := {[in,gn]}%=f denote the index sets identifying the line soli- 
tons as y — > 00 and as y — > — 00, respectively, according to Proposition ^. 51 Then two (A^_,A^ + )-soliton solutions of KPII 
are said to be in the same equivalence class if their asymptotic line-solitons are labeled by identical sets S± of index 
pairs, where \S+\ := N + =Nand |S_| :=A^_ =M-N. 

The set S+ US- of unique index pairings in Definition 12.71 has a combinatorial interpretation. Let [M] : = 
{1,2, . . . ,M} be the integer set with the pivot and non-pivot indices {ei, . . . ,ew} U {gi, . . . ,gM-N} forming a disjoint 
partition of [M]. Define the pairing map K : [M] — > [M] according to Proposition ^. 5f i) & (ii) as 

(2.7) %(e n ) =jn,n = 1,2,..., N, %(g n ) = i H , n = 1,2,...,M-JV. 

Then 7t : [M] — > [M] is a bijection, i.e., Jt G 5m, me permutation group of [M] [6|. In addition, Proposition 12.51 implies 
the following. 

PROPOSITION 2.8. The pairing map % defined by Eq. (|2.7| > is a derangement of [M] with N excedances, which 
are given by the pivot indices {e\, . . . ,eA?} of the coefficient matrix A in RREF. 

Recall that a permutation Jt with no fixed point is called a derangement, and an element / 6 [M] is called an 
excedance of Jt if Jt(/) > /. 

Each equivalence class of (N-,N+) -soli ton solutions of KPII is uniquely determined by a derangement Jt, as in 
Proposition ^. 81 These derangements also give a unique parametrization of a TNN Grassmann cell (see Remark 2.3) 
whose associated TNN matrix A satisfies the irreducibility Condition l2.4l Recall that in Example l2.6l both sets of line 
solitons as y — > ±°° are given by [1,4], [2,3], where 1,2 are the pivot indices and 3,4 are the non-pivot indices of the 
associated coefficient matrix A. In this case, the pairing map is a derangement of the set [4], and is given by 

(\ 2 3 4\ 
%= {4 3 2 lJ G54 ' 

in the bi-word notation of permutations in Sa- Note that the excedance set of Jt is {1,2}. In addition, Jt is also an 
involution of 54, i.e., Jt~' = Jt. Since the set of all involutions of Sin is isomorphic to the set of perfect matchings Mn 
introduced in Sectional, the involutions can be also represented by the chord diagrams representing the elements of 
SMjy. In particular, the chord diagram for the involution Jt e Sa given above depicts a nesting of the chords pi = [1,4] 
and P2 = [2,3] as shown below in Fig. 13.11 We remark that it is possible to represent derangements that are not 
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involutions by linear chord diagrams, with directed chords both above and below the line. These diagrams have been 
used to study the more general (7V_,7V + )-soliton solutions of KPII in Ref. [6|. But here we focus our attention to the 
7V-soliton solutions, which will be our next topic of discussion. 

3. 7V-soliton solutions 

When M — 27V, it follows from Proposition 12.51 that 7V_ = 7V+ = TV. If in addition, we consider S- = S+ in 
Definition 12. 71 then we recover the interesting subclass of the (N,N) -soli ton solutions mentioned in SectionQ] called 
TV-soliton solutions, which are characterized by identical sets of asymptotic line-solitons as |y| — > °°. Then the main 
features of the TV-soliton solutions follow primarily from our discussion in Section|2] in particular from Propositions |2.5l 
and !2.8l These are listed below. 

Property 3.1 TV-soliton solutions have the following properties. 

(i) The x-function of an TV-soliton solution is expressed in terms of 2TV distinct phase parameters and an TV x 27V 
coefficient matrix A which satisfies Conditions 12. 2| and |2.4| In addition, the TV x TV minors of A satisfy the 
duality conditions 1 6 , [201 : 

A(mi,...,m N ) =0 A(h,...,l N ) =0, 

where the indices {mi,... ,m^} and . . . ,lfj} form a disjoint partition of integers {1,2, . . .,27V}. That is, 
the phase combination 0(mi , . . . ,m#) is present in the T-function of Eq. d2.5t if and only if 8(/] , , Zjy) is. 

(ii) Each TV-soliton solution has exactly TV asymptotic line solitons as y — > ±°° identified by the same index pairs 
[en,gn] with e„ < g„,n = l,...,N. The sets {e\, ...,e^} and {g\ ,...,##} label respectively the pivot and 
non-pivot columns of the coefficient matrix A. Hence, they form a disjoint partition of the integer set [27V]. 

(iii) The amplitude and direction parameters of the n th asymptotic line soli ton [e„,g„] are the same as y — -> ±°°, 
and are given in terms of the phase parameters as 

fl « = kg n — k en , c„ = k gn + k £n . 

(iv) The pairing map associated to an TV-soliton solution, namely n(e„) = g n , n(g n ) = e„, n= 1,2, ... ,7V, corre- 
sponds to a partition of the integer set [27V] into TV distinct pairs of integers (e n ,g„), as in Sectional- Each 
such map is an involution in Sin with no fixed points, a member of 

i2N = {n e Sin : = JtandTi(z') ^ ;', V/ 6 [27V]} . 

Such permutations can be expressed as products of TV disjoint 2-cycles, and their chord diagrams are iden- 
tical to those of the perfect matchings Mn (see Fig. 12. U . The total number of such involutions is given by 
\Iw\ = \Mn\ = (2TV — 1)!!. Hence, there are (27V — 1)!! distinct equivalence classes of TV-soliton solutions. 

3.1. Equivalence classes of 2-soliton solutions. When TV = 2, there are three types of 2-soliton solutions referred 
to as the O-, T- and P-types (following the terminology introduced in Ref. [20|). They are identified by the canonical 
coefficient matrices associated with T-functions, namely 

(\ 1 0\ , (\ -1 -A , (\ -A 

Ao = {o o i ij' ^={o i Xl X2 )> Ap= {o 1 1 oj' 

with xi > X2 > in At- By applying the rank conditions of Proposition ^. 5l to the above coefficient matrices, it is easily 
verified that the O-type 2-solitons have asymptotic line-solitons [1,2] and [3,4]; the T-type resonant 2-solitons have 
asymptotic line-solitons [1,3] and [2,4]; and the P-type 2-solitons have asymptotic line-solitons [1,4] and [2,3]. These 
are shown in Fig. 13. II Notice that each of the O- and P-type solitons interact via an X-junction. After interaction, each 
line soliton undergoes a position shift in the xy-plane. However, it can be shown that the position shifts for the O-type 
solitons are opposite in sign to that of the P-type solitons (6l . On the other hand, the T-type solitons interact via four 
Y-junctions, connecting the four asymptotic line-solitons to four intermediate segments. Each of these intermediate 
segments is also a line soliton. For example, in the T-type solitons in Fig. 13. li the asymptotic line soliton [1,3] (as 
y — > — oo) forms the intermediate line-solitons [1,2] and [2,3] at the bottom left Y-junction. The line-soliton [2,3] 
connects with the asymptotic line-soliton [2,4] (as y — > <*>) and the line-soliton [1,2] connects with the asymptotic 
line soliton [2,4] (as y — > — oo). Similarly, the asymptotic line soliton [1,3] (as y — > oo) forms the intermediate line- 
solitons [1,4] and [3,4] at the top right Y-junction. The line-soliton [3,4] connects with the asymptotic line-soliton 
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O-type 



T-type 



P-type 



Figure 3.1: Three different two-soliton solutions of KPII with the same phase parameters (fej , 
the three 2-soliton equivalence classes: O-type, T-type and P-type 2-soliton solutions. 



fct) = (—2,— A, 0,1), illustrating 



[2,4] (as y — » oo) and the line-soliton [1,4] connects with the asymptotic line soliton [2,4] (as y — > -«>). Fig. 13.11 
also shows the chord diagrams for the corresponding pairing maps, which are involutions of the permutation group 
54. They correspond to the disjoint partitions of [4] into 2 pairs. In cycle notation, these involutions are given by 
7io = (12) (34), jct = (13) (24), and ftp = (14) (23), for the O-, T- and P-type 2-soliton equivalence classes. According 
to the Definition l2.il the chord diagram for Tto forms an alignment; whereas the diagram for jit has a crossing between 
the chords corresponding to the line solitons [1 , 3] and [2,4]; and the diagram for Tip is a nesting. 

An important distinction among the three types of 2-soliton solutions is that they belong to different regions of 
the soliton parameter space. Suppose (ai,ci) and {a%,C2) are the soliton parameters of the asymptotic line-solitons of 
each type, with the same set of distinct phase parameters. Since the phase parameters are ordered: k\ < ■ ■ ■ < k\, the 
soliton parameters satisfy the following relations, which can be easily verified using Eqs. d3.lt . 

(i) For O-type 2-soliton solutions, C2 > c\ and C2 — c\ > a\+ai- 

(ii) For T-type 2-soliton solutions, cj_ > c\, and \a\ — ci2\ <C2 — c\ < a\+ ai- 

(iii) For P-type 2-soliton solutions, > fli and \c2 —c\\ < — a.\. 

(iv) (c2-ci)o > (c2-ci)t > |c2-ci|p, (oi + 02)0 < («i + ai)T = (fli + ai)v , and 
|fl2-ai|o = ks-aik < (fl2-ai)p- 

Note that for O- and T-type solutions the soliton directions are ordered, while for P-type solutions the amplitudes 
are ordered. Any choice of the soliton parameters a\,c\\a%,C2 with 01,02 > would lead to one of the three types 
of 2-soliton solutions, provided that {c\ ± a\ , C2 ± 02} are distinct real numbers. Thus, the three types of 2-soliton 
solutions partition the soliton parameter space into disjoint sectors, bounded by the hyperplanes \c2 —c\ \ =01+02 
and \c2 — c\ \ = \a\ — 02!- At each boundary between two sectors, two of the phase parameters coincide. In such a 
situation, it can be shown (by taking suitable limits) that the 2-soliton solution degenerates into a Y-junction Il5l l20l . 

It should be clear from the above that the O-, T- and P-type 2-soliton solutions exhibit distinct types of interaction 
patterns, and belong to different regions of the soliton parameter space. For > 2, in addition to the non-resonant (O- 
and P-type) and fully resonant (T-type) solutions, a large family of partially resonant solutions exists. For example, 
when N = 3, Property l3.ir iv) implies that there are 15 distinct equivalent classes of 3-soliton solutions (see Fig J3.2b . 
Unlike the N = 2 case above, it turns out to be a complicated task to classify the A^-soliton solutions according to 
their coefficient matrices A, when N > 2. This task was recently carried out by the authors, and will be reported in a 
future publication J7|. Here we consider a more direct classification scheme for the A^-soliton equivalence classes, by 
characterizing the pairwise interactions between the line solitons of O-, T- and P-type, much as in the 2-soliton case. 
In other words, we represent the A^-soliton solutions by the corresponding involutions in I2N C S2N (equivalently, the 
matchings in Mm), and enumerate the solutions according to the number of alignments, crossings and nestings of the 
associated chord diagram. We describe this classification scheme below. 



3.2. Combinatorics of A^-soliton solutions. Throughout this subsection, we associate an A^-soliton equivalence 
class defined by the set S — {[e„,gn}}^ = i of asymptotic line solitons (see Definition 12.71 ) with the partition X — 
{pi,P2, ■ ■ ■ ,Pn} £ Mn, where p„ :— [e n ,g n ]. Recall from Property |3.1f ii) that the integer set [2N] is a disjoint union of 
the index sets E := {e\, . . . ,ew} and G := {gi, . . . with the following orderings among the indices: 

(i) 1 = e\ < e 2 < ■ ■ ■ < e N < 2N, 

(ii) e n < g n for all n — 1,2, . . . ,N. 

An immediate consequence of the above orderings is that 

(3.2) n < e n < 2n — 1 , n — \,...,N, 

since there are at least n — 1 indices to the left of e„, namely ei,e2, • • • ,e n -i', and at least 2N — 2n + 1 indices to the right 
of e n , namely g n ,e r ,g r , r > n. The A^-soliton classification scheme is obtained by considering various statistics over 
the possible chord configurations for the chord diagrams of Mn- For this purpose, using Definition 12. II we introduce 
the following sets, which record the total number of alignments, crossings and nestings for a given chord in any chord 
diagram of Mn- 

Definition 3.2. Let p„ = [e„,g„] be a given chord of a partition X € 5Wa?, and let B n := {p r = [e r ,g r ] : r <n} 
be the subset of chords originating from the left of p„ in the linear chord diagram of X. 

(a) The set O n of alignments with the chord p„ forming O-type configurations, and the alignment number 
al(X), are defined by n 

O n := {p r = [e n gr] S B n : g r < e n } , al(X) := ^ \0„\. 

n=\ 

(b) The set T„ of crossings with the chord p„ forming T-type configurations, and the crossing number cr(X), 
are defined by n 

T„ := {p r = [e r ,gr] 6 B n : e„ < g r < g n ] , cr(X) := ^ |7;,| . 

n=l , 

(c) The set P„ of nestings with the chord p„ forming P-type configurations, and the nesting number ne(X), are 

defined by r , * 

Pn ■= {Pr = [e r ,g r ] G B„ : g r > gn) , ne(X) := Y \Pn\- 

11=1 

It follows from the above definitions that B„ is the disjoint union of the sets O n , T n and P„, so that \O n \ + \T„\ + 
\P„ \ =n—l, and al(X) +cr(X) +ne(X) = A^(A^ — l)/2, which is a count of all possible pairwise chord configurations 
in the partition^. Note that for 0„, the indices g r lie in the intervals (e r ,e r+ i), 1 < r < n. Hence, \0„\ = e„ —n. So 
the number of crossings and nestings with the chord p„ sum to \T„ \ + \P n \ = (n — 1) — (e„ —n)=2n — e n — \, which 
depends only on the pivot index e n € E. This observation leads to the following. 

LEMMA 3.3. If fM (E) C denotes the set of all partitions which have the same (pivot) index set E, then the 
number of partitions of 94 (E) having r crossings and s nestings is the coefficient of p s q r in 

N p" —q" 

m E (p,q) = Y[ [2n-e„] P 4, [n] p . q := — - — = £ p'q 1 . 

n=\ P 1 i+j=n-\ 

The degree of both p and q in mg (p, q) is N 2 — {e \ + e 2 + . . . + e^). 

Proof. The distribution of crossings and nestings is the sum of p ne ( x ) q cr ( x *> over all partitions X G M (E). Using 
Definition l3.2l for cr(X) and ne{X), this distribution can be expressed as 

N 2n—e n —\ 

Y p(\ p W\+-+\P(N)\) q (\TW\+-+\T(N)\) = ^ Y p l q 2n ~ e "~ l ~ l , 
Xe»t(E) n=l 1=0 

after interchanging the sum and product, and using the fact that \T n \ + \P n \ = 2n — e„ — 1 for n = 1 , 2, . . . , N. Since the 
second sum is precisely [2n — e„} p q , the formula for m^(jp^q) follows. □ 

It is easy to verify from the product formula that mE(p,q) is symmetric in p and q. Consequently, the number of 
diagrams with r crossings and s nestings is the same as the the number of diagrams with s crossings and r nestings 1 19 1. 
Note also that the enumerating polynomial for the crossings alone is given by m^Lg); while niE(p, 1) enumerates 
only the nestings for the chord diagrams of M (E). In order to extend the results of Lemma [3~3l to the entire set 
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Mn, one needs to sum niE(p,q) over all possible choices of the integer set E, with e„ E E satisfying Eq. (13.2b . Using 
Lemma [3~3l the expression for the required generating polynomial is given by 

(3.3) F N (p,q):= £ p^V^ = J> £ (p, 9 ) = E f[[2n-e n ] p>q . 

XeM N {£} l=e l <e 2 <...<e N , n=l 

k<e k <2k-\ 

Since we have from Eq. d3.2t that e n > n, n = 1,2, ...,N, it follows from Lemma [331 that the degree of p and g in 
F N (p,q) is given by 

2 „ ■ „ , , .n.^- 1 ) 



W = cr(X)m« = iV z - (1 + 2 + . . . + N) 
Furthermore, like rriE(p,q), Fff(p,q) is symmetric in p and q, i.e., 

N(N-l)/2 

F N (p,q) = £ c rs q r p s , c rs = c sr . 

r,s=0 

Some interesting consequences of Eq. (13. 3t for special cases of F^(p,q) are collected below. 
COROLLARY 3.4. The function FN(p,q) has the following properties: 

(i) 2^(1,1) = |^v| = (2JV-l)H. 

(ii) When p = 1 , q = 0, f/ie tofa/ number of non-crossing ( i.e., only alignments and nestings) chord diagrams of 
!Mn equals the N th Catalan number [ 10 1. That is, F ( 1 , 0) = Cn which also counts the total possible choices 
for the ordered integer set E [6 |. Similarly, Fn(0, 1) = Cn gives the total number of non-nesting (i.e., only 
alignments and crossings) chord diagrams of Mm- 

(iii) Fff(l,q) —: Fff(q) gives the generating polynomial for the number of crossings introduced in Section\2\l , 
which is given explicitly by the Touchard-Riordan formula Eq. ( 12. lb . 

The polynomials Fp{(p,q) can be determined from a generating function F(p,q,x) which is a formal power series, 
and has the following representation. 

PROPOSITION 3.5. The generating function for FN(p,q) is the Stieltjes-type continued fraction, namely 
F(p,q,x) := £ F N (p,q)x N = — , F (p,q) := 1 . 

N=0 j x i>-\p,q 



x[2] p , c 



1 — 

Proof. First consider the set E := {1 = e\ < ■ ■ ■ < e^ : e^ < 2k — 1}. Note that E can be decomposed into distinct 
subsets when e„ =2n—l. One has 

2V-1 

E= U (E n UE n ) , 

n=0 

where E n :={\ = e\ < ■ ■ ■ < e„ : e^ < 2k — 1 } for n ^ can be viewed as the n-truncates of the original set E,Eq = 0, 
and 

E n := {2n + 1 = e n+ \ < ■ ■ ■ < e N : 2n + k < e„ +k < 2(n + k) - 1} . 
The set E„ can be re-expressed as 

E' N _ n = {l = e\<-~<e l N _ n :k<e' k <2k-l}, 

by shifting and relabeling the indices as e n+ k : = e' k + 2n. Note however that E' N _ n (with all e' k < 2k — 1 ) is not the same 
as £W_„ (with all e k < 2k — 1). 

From Eq. (13.31 >. 

N N-l 

(3.4) F N (p,q) = £ J! \ 2k - e k] P ,q = £ F n(P,l) c N-n(p,q), 

{E} k=\ n=0 
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where C n (p,q) = £ [2k — e^] M . Introduce the power series 



F (PAi x ) = £ f n(pa)^ , F (p,q):=l. 



N=0 

C(p,q,x) = £ C N {p,q)x N . 
N=\ 

Using Eq. ( 13.41 ) in the power series, one finds that F(p,q,x) — 1 equals the product F(p,q,x)C(p,q,x), which implies 

(3-5) F(p,q,x) = — I r . 

1 -C(p,q,x) 

Next, define the associated polynomials 

n n 

F n (p,q;l) := £ Y\[ 2k ~ e k + l] P ,q , C n (p,q;l) := £ Y\[2k-e' k + l] p , q , 

{E n }k=l {E'}k=l 

so that F„(p,q;Q) = F n (p,q) and C n (p,q;0) = C„(p,q). The corresponding power series F(p,q;l,x) and C(p,q;l,x) 
are defined similarly to F(p,q,x) and C(p,q,x) above, and they also satisfy Eq. d3.51 l. Furthermore, for n > 1 the 
associated polynomials satisfy the relation 

(3-6) C„( P ,q;l)= £ Il[2A:-4 + 'W=[' + l] M E fl[2£-4 + /U 

{E>}k=l {E' n }k=2 
n-\ 

(3.7) =[/+iu I n[ 2 ^'-^+( / + 1 )W = [ /+1 U F »-iO'^; / + 1 )' 

after an appropriate index shift, k = j+l, and relabelings e'- + : = e ; - + 1 so that j <ej <2j — 1 for /' = 1 , . . . n — 1 . As 
a result, the set £^ changed to the set E n -\. The formal power series constructed from the first and last expressions in 
Eq. (13.7b satisfies C(p,q;l,x) — x[l + \ ] p q F(p : q;l + l,x). From the analogue of Eq. (13.5b for the associated functions 
F(p,q;l,x) andC(p,q;l,x), one therefore obtains 

1 

F(p,q;l,x) 



1 - [l+l] P;q xF(p,q;l + l,x) 
This yields the continued fraction representation for F(p,q,x) = F(p,q;Q,x). □ 



We can graphically illustrate the results of Proposition [33] for = 2 and 3 in terms of the corresponding chord 
diagrams. For N = 2, F2(p,q) = 1 + p + q, which implies that there is one each of the O-, T-, and P-type diagrams. 
These were displayed in Fig. 13.11 For N ~ 3, there are 15 chord diagrams, which are displayed in Fig l3.2l They are 
characterized by 

F 3 (p,q) = (l+2p + p 2 + p 3 ) + (2 + 2p + 2p 2 )q+{l+2p)q 2 + q 3 . 

Note that the total number of pairwise chord configurations for each case for = 3 is 3(3-l)/2 = 3. We use the ordering 
(P1P2 1 P2P3 1 P3P1 ) for the chord-pairs, and indicate the interaction type for each pair below: 

(a) Non-resonant cases: 1 (PPP)-type, 1 (PPO)-type, 2 (POO)-type, 1 (OOO)-type. These are on the first row 
in Fig l3.21 from left to right. 

(b) One-resonant cases: 2 (TPP)-type, 2 (TPO)-type, and 2 (TOO)-type. These are on the second row. 

(c) Two-resonant cases: 2 (TTP)-type and 1 (TTO)-type. These are on the the third row. 

(d) Three- (i.e., fully-) resonant case: 1 (TTT)-type. 
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(POO) 



(OOO) 



(TIP) 




(TOO) 



(TTP) 



(TTO) 




(TTT) 



Figure 3.2: The closed chord diagrams for 3-soliton solutions. The dots indicate the pivots (ei,e2i e 3), an d the ordered letters 
below each diagram indicate the type of interactions in (piP2,P2P3,P3Pl) with the soliton pairing p„ = [e n ,g n ]. The number of the 



diagrams having the same number of crossings comes from the generating function Fj(q) 
number C3 = ^3(0), counts the diagrams in the first row. 



■ q } + 3q l + bq + 5. E.g., 5, the Catalan 



3.3. Generating function and ^-orthogonal polynomials. In the special case when p = 1, the formula for 
F(q,x) := F(l,q,x) in Proposition 13.51 reduces to similar continued fraction expression for F(q,x) := F(l,q,x), 
namely 



(3.8) 



F(q,x) = £ F N (q)x N 

N=0 



1 



--[2]. 



:[3], 



Fo(q):=h 



with [n] q := 1 + q + ... + q"^ 1 . From Corollary I3.4t iii). it follows that F(q,x) is the generating function for the 
polynomials Fm(q) that enumerate the crossings of the chords of Mn, whose explicit formula is given by Eq. ( 12. jj . 
Here we show that the continued fraction for F(q,x) is related to the moment generating function for the continuous 
g-Hermite polynomials. It follows that the Touchard-Riordan polynomials F^iq) are simply the even moments of the 
weight function with respect to which the g-Hermite polynomials are orthogonal. The latter result was also found in 
Ref. 03. 

We first collect some facts (see e.g. |HJ[21]) from the spectral theory of bounded, real, semi-infinite Jacobi matrices 
on the Hilbert space 



fe(C) := {u = (mq, wi ,U2, . . .) : u, 6 ( 
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k=0 



< 00} . 



Define the following tri-diagonal matrices 



/0 1 

a\ 




1 



a 2 



with a,- > 0, i = 1,2, . 



fO 1 

fli 





1 







/0 1 

02 



•■■ fl„_ 2 1 ■•■ 

Next, consider the linear system of equations 

(Xl-L)(f> = e , for 9= (90, 91,92,... ) r , 

where / := diag(l, 1, . . .) is the semi-infinite identity matrix and eo = (1,0,0, . . .) £ /2(C). 

Fact (a) (Resolvent of Jacobi matrix): The (0,0)-element of the resolvent of L, i.e., 90 = 
developed in a continued fraction by rewriting the linear system of equations as follows: 

1 

90 = and 




1 



a n -2 








1 

a n -\ 0/ 



(XI — L) eo), can be 



9» 



n > 1 . 



9o 



9«-i 



X- 



9»+i 
9« 



Thus one has 
(3.9) 



9o : 



X-- 



«2 



Forn > 0, then convergent of this continued fraction is of the form^„ = — , where No = 0, Do = 1 andA^„(A,), D n (X) n > 

1, are polynomials in X of degrees n — 1 and n, respectively. These sequences of polynomials satisfy the 3-term recur- 
rence relation 

(3.10) XP„ = a n P n - y+Pn+u P n := (N n , D n ) , 

with Pq = (0, 1) and Pi = (l,X) as respective initial conditions. Furthermore, it can be shown by Cramer's rule that 
for n > 2, 

D„ = det(Xln-L n ), Nn = 6etQJ n -i—L n -i), suchthat ^ = (e , (XI „ -L„) _1 e ) , 

L) n 

where /„ is the nxn identity matrix and eL = (1 , . . . , 0) G K". 

Fact (b) (Spectral theorem): If {a n }, n > 1, is a positive bounded sequence such that the Jacobi matrix L, defined 
above, is bounded on /2(C), then there exists a unique spectral measure /j with compact support E such that 



(3.11) 



90 = lim ^ = (e 0) (XI — L)~ 1 eo) = [ ^ , for X £ E 



d/j(s) 



Furthermore, the polynomials {N n , D„},n>l, are orthogonal with respect to the measure /j. In particular {£>„}, n > 1 
form a sequence of monic polynomials satisfying the orthogonality relations 



/ D m (s)D n (s)dp.(s) = a„8 m „ , with a„ = J~[a ; - 



,/=i 

If we now set a n — [n] q for n > 1 in the continued fraction representation for 90 in Eq. ( 13.91 ), and compare the 
resulting expression with the generating function in Eq. (13.81 ). we find that 

F(q,x = X )=^— - 

Af=0 ^ jfc=0' 

where the last equality follows from Eq. ( 13.111 ). Note also from Eq. ( 13.111 ) that the moments J L s k dp.(s) = (eo, L k eo) , k = 
0,1,2,..., clearly vanish for any odd k because of the structure of the Jacobi matrix L. Therefore, we conclude that the 

13 



generating polynomials Fp](q) are given by the even moments of the measure p.. Moreover, Eq. ( 13.101 ) with a n = [n] q 
is the well known 3-term recurrence relation for the g-Hermite polynomials H n (s, q) (see e.g., Hm ). Indeed, it follows 
from the initials conditions that the denominator polynomials above satisfy D n (s) = H„(s,q) for ra = 0, 1,2, .... These 
polynomials satisfy the orthogonality relations 



J — i 



H m (s,q)H„(s,q)dp(s) = [n] q \h mni 

dp(s,q) :=v(s,q)ds , sE[—a,a], a : = 2(1 — q)~ l l 2 , 
where \q\ < 1, and the weight function v(s,q) is given in terms of the Jacobi theta function ©i (Q,q) by 

_ 1 ^ 

V(*,g) - ^-®i(B,q) = 2- £(-l)V ( " +1)/2 sin(2» + l)9, 

cos9:=- = ±(l-4) 1/2 s. 
a 

Note that v(s,q) is an even function in s (i.e., it is stable under 8 — > Jt — 8, where 8 G [—71, jc]), so the odd moments 
vanish as observed above. The even moments for the g-Hermite polynomials are given by 

/a 
x 2N v(q,x)dx 
-a 

r>2N 



£* ( - 1 f q < n+1 )l 2 r cos 2N 8(cos 2«8 - cos2(n + 1 )6) dd , 
,_n JO 



which yields the Touchard-Riordan formula, after evaluating the last integral and rearranging the summation indices. 

Remark 3.6 It is intriguing to note that the generating function F(p,q,x) that enumerates the interaction types of the 
Af-soliton solution of the KPII equation has its origin in the theory of ^-orthogonal polynomials. We mention another 
relation with orthogonal polynomials, without presenting the details. Consider the case p = 0. It turns out that the 
generating function F(Q,q,x) for the non-nesting chord diagrams is related to the moment generating function for a 
certain class of ^-orthogonal polynomials studied in Ref. [2| (see also [16]). A particularly interesting consequence of 
this relation is that the function F(0,q, —q) has a Rogers-Ramanujan interpretation. Let (pi and 92 be certain modular 
forms of weight I for the level-5 principal modular group T(5) < PSL(2,Z); namely, 

(3.12) = * £(-1)^0^/40, ^ = 1 

where r\(q) — q 1 ^ 24 ]~[(1 — q") is the Dedekind r|-function. It is well-known that the modular forms (pi and 92 admit 

n=l 

infinite product representations, which constitute the Rogers-Ramunajan identities. Accordingly, F(Q,q, —q) can be 
represented as a quotient: 

1/5 92 _ fr tl-ntl-^) _ 1 



F { 0,q,-q) = q-y^ r m 



9i ,i A o(l-9 5 " +2 )(l-? 5 " +3 ) 



? 2 



? 3 



This is the Rogers-Ramanujan continued fraction 1 27l l29l . 



4. Conclusion 



We have presented a classification scheme for the A^-soliton solutions of the KPII equation, based on the combi- 
natorics of chord diagrams consisting of chords connecting distinct pairs of 2N points. We have shown that it is 
possible to associate the O-, T-, and P-type of pairwise interaction patterns among the asymptotic line solitons of the 
Af-soliton configuration with the alignments, crossings and nestings among pairs of chords in the chord diagram. As a 
result, the equivalence classes of A^-soliton solutions can be enumerated by the same generating polynomial Fn{p,q) 
(Eq. (13.3b ) of the distribution of nestings and crossings for the set Mn of all chord diagrams of [2N]. It follows from 
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Propositions 2.5 and 2.8 that each asymptotic line soliton of a given A/-soliton solution is uniquely identified with an 
index pair p„ = [e„,g„] , 1 < e„ < g„ < 2N, which also labels a particular chord in a chord diagram. This pairing map 
(Eq. j2.1V ) plays a crucial rule in establishing a correspondence between an AZ-soliton equivalence class and a particu- 
lar chord diagram of 9vt^. The soliton pairing encoded in Proposition 2.5 can be derived via a systematic asymptotic 
analysis (for fixed t ) of the A/-soliton T-function, which is sum of real exponentials with positive coefficients, by identi- 
fying those phase combinations ®{m\ , . . . , wijy) that are dominant in different regions of the xy-plane as \y \ — > °° Since 
a discussion of the asymptotics of the T-function is beyond the main focus of the present article, it has been omitted 
here. Interested readers may find the relevant details, including a proof of Proposition 2.5, in (3] (see also 151). 

Finally, we have derived a continued fraction representation for the generating function F(p,q,x) of the poly- 
nomials FN(p,q). We have shown that the special cases of this generating function are, in fact, moment generating 
functions of certain kinds of g-orthogonal polynomials. It is interesting to speculate whether the unrestricted AZ-soliton 
generating function F(p,q,x) is the moment generating function for some new family of (p, g)-orthogonal polynomi- 
als. 
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